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Abstract 

We study if leptogenesis works successfully together with the neutrino mass anarchy hypothe- 
sis. We find that the predicted neutrino mass spectrum is sensitive to the reheating temperature 
or the inflaton mass, while the distributions of the neutrino mixing angles and CP violation 
phases remain intact as determined by the invariant Haar measure of U(3). In the case of ther- 
mal leptogenesis, the light neutrino mass distribution agrees well with the observations if the 
reheating temperature is O{10^~^^) GeV. The mass spectrum of the right-handed neutrinos and 
the neutrino Yukawa matrix exhibit a certain pattern, as a result of the competition between 
random matrices with elements of order unity and the wash-out effect. Non-thermal leptogen- 
esis is consistent with observation only if the inflaton mass is larger than or comparable to the 
typical right-handed neutrino mass scale. Cosmological implications are discussed in connection 
with the 125 GeV Higgs boson mass. 
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I. INTRODUCTION 



The origin of flavor of quarks and leptons in the standard model (SM) remains one 
of the great mysteries in the particle physics: why are there three generations of quarks 
and leptons, and why hierarchical structure in the Yukawa couplings? There have been 
proposed a variety of models, some of which rely on a hypothetical flavor symmetry. 
However, none of them are decisive as yet. This is partly because the mass spectrum of 
elementary particles shows no definite pattern, unlike the periodic table of elements, and 
so, if we are to understand the flavor structure based on the fundamental symmetry prin- 
ciple, either small symmetry breaking or unknown coupling constants must be introduced 
ad hoc, allowing a great variety of flavor symmetries and charge assignments. 

While symmetry has been a useful and attractive guiding principle in physics, it is 
not necessarily applicable to all observables. For instance, the observed vanishingly small 
cosmological constant may be interpreted to suggest the existence of some profound fun- 
damental symmetry setting the cosmological constant (almost) zero. Alternatively, it may 
simply be that the cosmological constant is an environmental parameter adjusted by the 
anthropic principle Similarly, some of the observables in our Universe may be strongly 
affected by the anthropic conditions, and if so, it is hopeless to try to understand their 
values from the symmetry principle. Indeed, many parameters in the SM seem to be 
adjusted so that the existence of life is possible. Therefore, the apparent pattern of the 
quark and charged lepton mass spectrum may be just a consequence of the environmental 
selection, and may not reflect any fundamental symmetry.^ 

On the other hand, the situation is different in the neutrino sector. The neutrinos are 
massless in the SM, but the neutrino oscillation experiments have revealed that neutrinos 
have a tiny but non-zero mass j^. Its typical mass is constrained to be below 0.2 eV jsj, 
and its cosmic mass density is much smaller than the observed dark matter density. The 

^ That said, it is difficult to experimentally confirm such anthropic argument. In the rest of this paper 
therefore we do not attempt to interpret the structure of Yukawa couplings for quarks and charge 
leptons. We will come back to the flavor symmetry later. 
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tiny neutrino mass can be beautifully explained by the celebrated see-saw mechanism ^; 
the smallness of the neutrino masses is related to the ratio of the weak scale and the heavy 
right-handed neutrino mass Mq ~ 10^^ GeV close to the GUT scale. With such tiny mass 
and cosmic energy density, therefore, the neutrino mass and mixing may be irrelevant to 
the existence of life, and so, it may possess information on its original distribution in the 
landscape. 

Let us briefly summarize the current status of the neutrino parameters. The three 
neutrino mixing angles are given by jsj: 

sin2 0i2 = 0.32010:°J57, 

sin^e^s = 0.49l°:°« (0.53lH^), 

sin^^i3 = 0.026^^04 (0.027^°:°°^), (1) 

where the normal (inverted) hierarchy is assumed. We note that two of them are large, 
but even the smallest one, ^13, is not extremely small. The mass squared differences are |5| 

Am^i = (7.62 ± 0.19) X 10~^eV^ 

Ami, = 2.53^Ho i-^-'^^-om) x 10"' eV^. (2) 

The ratio of the mass squared differences is Am2i/|Amg|| ^ 0.03, which is much milder 
compared to that for the quarks or charged leptons [6]. Intriguingly, those neutrino 



parameters are consistent with the neutrino mass anarchy hypothesis proposed in Ref 
which has been further studied in Refs 



7| 



4lOl|. In particular, the recent discovery of 
non-zero ^13 by Daya-Bay experiment ll| has made the idea very attractive [lO^.^ As we 
shall briefly review in the next section, the neutrino mass anarchy hypothesis is based on 
two assumptions: (i) there is no quantum number to distinguish flavors in the neutrino 
sector, and therefore the couplings are structureless in the flavor space; (ii) the couplings 
and mass matrix obey basis-independent random distribution. In particular, it was shown 
in Ref. [sj] that the mixing angle distribution obeys the invariant Haar measure of U(3). 
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A hint for non-zero 6*13 was reported by T2K 12 1, MINOS 11311 and Double-Chooz [14!] experiments. 



Recently, RENO experiment also observed the non-zero 613 [15[. 
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Also, using the linear measure of the eigenvalues of the random matrices, the observed 
neutrino mass squared difference can be naturally explained in the see-saw mechanism. 
Thus, both neutrino mass anarchy and the see-saw mechanism are arguably the most 
attractive framework for understanding the observed neutrino parameters. 

The origin of matter remains a puzzle in cosmology and particle physics. Since any pre- 
existing baryon asymmetry would be exponentially diluted by the subsequent inflationary 
expansion, it is necessary to generate the baryon asymmetry after inflation. One plausible 
explanation is the baryogenesis through leptogenesis 16| : the lepton asymmetry generated 
by the out-of-equilibrium decay of the right-handed neutrinos is transferred to the baryon 
asymmetry via the sphaleron process. However, if the leptogenesis is responsible for the 
observed matter asymmetry, it might select a certain subset of the neutrino parameters, 
and as a result, the original distribution in the landscape may be significantly distorted, 
spoiling the success of the neutrino mass anarchy. 

In this paper we study if the leptogenesis works successfully together with the neutrino 
mass anarchy hypothesis. The result is two-fold. First, the mixing angles and the CP vio- 
lation phases (one Dirac and two Majorana) in the low energy are subject to the invariant 
Haar measure of U(3), and they are not modified by requiring the successful leptogenesis. 
In a sense, the mixing angles as well as the CP violation phases are orthogonal to lep- 
togenesis. Second, the neutrino mass eigenvalues are generically affected by leptogenesis. 
We find however that thermal leptogenesis is possible without significant modification 
of the predictions of the original neutrino mass anarchy, if the reheating temperature is 
O{10^'^^) GeV and if the typical right-handed neutrino mass scale is of O{10^^) GeV. This 
is the result of the competition between random matrices of order unity and the wash-out 
effect. As a result, the mass spectrum for the right-handed neutrinos and the neutrino 
Yukawa matrix exhibit a certain pattern, which is quite similar to that can be understood 
in terms of a conventional flavor symmetry. In other words, the flavor symmetry of the 
right-handed neutrino sector is emergent in this framework. In the case of non-thermal 
leptogenesis, we flnd that the neutrino mass spectrum is signiflcantly affected in contra- 
diction with the observations, if the inflaton mass is smaller than the typical right-handed 
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neutrino mass scale of 10^^ GeV. It suggests that the inflaton mass needs to be larger 
than or comparable to O{10^^) GeV for successful non-thermal leptogenesis. 

The rest of this paper is organized as follows. In Sec. [TTlwe briefly review the neutrino 
mass anarchy hypothesis and define our notation and framework. In Sec. IIIII we discuss 
how leptogenesis affects the neutrino parameters. The last section is devoted for discussion 
and conclusions. 



II. NEUTRINO MASS ANARCHY 



In this section we briefly review the neutrino mass anarchy and its prediction, focusing 
on the see-saw mechanism. 



A. Preliminaries 



We consider the following see-saw Lagrangian, 

C D emijH + hij NdjH - ^M.jNiN, + h.c, (3) 

where i,j = {1,2,3} denote flavor indices, ii represents the left-handed lepton doublets, 
em are the charged lepton singlets, Ni are the right-handed neutrinos, and H is the Higgs 
doublet, fij and hij form complex- valued 3x3 matrices of charged- lepton and neutrino 
Yukawa couplings, respectively, and Mjj forms a complex valued 3x3 symmetric matrix 
of the right-handed neutrino Majorana mass. For later use we also define a dimensionless 
matrix Xij as 

X, ^ ^, (4) 

where Mq is the typical mass scale of the right-handed neutrinos. Mq can be interpreted 
as the B— L breaking scale, if the right-handed neutrino mass arises from the vacuum 
expectation value (vev) of the B— L Higgs boson through a renormalizable interaction 
with a coupling of order unity. We adopt Mq = 10^^ GeV as a reference value throughout 
this paper. Later we will briefly discuss how our results will change for different values of 
Mo. 
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Quantum mechanics dictates that any states with the identical quantum numbers 
should mix with each other. If there is no quantum number which distinguishes three 
generations of ^i and A^,, the matrices hij and Mjj are considered to be structureless. In 
particular, they may be subject to a basis-independent random distribution. This is the 
essence of the neutrino mass anarchy. On the other hand, the charged lepton mass matrix 
(as well as that for quarks) is probably determined by other physics such as the anthropic 
considerations or conventional flavor symmetries, and so, we do not attempt to interpret 
the structure of fij in terms of the anarchy here (see e.g. 17Hl9l|). Therefore we simply 
adopt a basis where the charged-lepton Yukawa matrix is diagonalized: 



H + h,MaH - -MijN^Nj + h.c. 



(5) 



where a and (3 represent the lepton flavor indices, e, /x, r. The anarchic nature of hia and 



Mij is maintained in this basis. 



The neutrino Yukawa matrix {h)ia can be diagonalized by the bi-unitary transforma- 
tion. 



h UlhUi 



^ h ^ 



V 



/l2 

h 



(6) 
(7) 

(8) 



where Ul and Ur are unitary matrices, and we take < hi < h2 < h^. Similarly, the 
right-handed neutrino Majorana mass matrix can be diagonalized as 



31 



^ Ml Q ^ 



v 



Ms 
M3 



= D 



(9) 
(10) 



where Ujsf is a unitary matrix, and one can take < Mi < M2 < M3 without loss of 
generality. 
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Below the scale of the right-handed neutrino, we obtain a low-energy effective interac- 
tions containing a Majorana mass for left-handed neutrinos: 



(11) 



where 



K)a/3 = [h^X-'h]^ 



' Mo 

= (uiDr^UlUl^D-^lU^UnDnUl) v\ (12) 

with i; ~ 174 GeV is the vev of the Higgs field. The light neutrino mass can be naturally 
explained by the heavy right-handed neutrino mass Mq ~ 10^^ GeV close to the GUT 
scale in the see-saw mechanism j4|. 

The neutrino mass can be diagonalized as 



mi ^ 



11* 

^ MNS 



\ 



u 



t 

MNS^ 



(13) 



ms 
mg 

where Umns is a unitary matrix. There is currently no constraint on the sign of Am'^^, but 
a neutrino mass spectrum with normal hierarchy is preferentially realized in the neutrino 
mass anarchy, and so, we will assume < mi < m2 < unless otherwise stated. 

Note however that, although rare, the inverted mass hierarchy is possible. But one 
should be careful when comparing the result with the observations ([T]) and ([2]), because 
ms is always the heaviest in our notation. As far as the mass difference squared is 
concerned, one should simply replace mi — )■ m^, m2 — > mi and — )■ m2 in order to 
compare our results for the inverted hierarchy with the observations. We will come back 
to this issue at the end of this section. 

The neutrino mixing matrix Umns can be parametrized as follows. 



U 



MNS 



C12C13 



S12C13 



-S12C23 - £125235136*"^ C12C23 - Sl2S23'Sl3e*^ 



S12S23 - Cl2C23Sl3e 



-C12S23 - Sl2C23Sl3e 



Sise"'" 
S23C13 
C23C13 



X diag 1, e 2 ,62 



(14) 
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where = cos^, Sij = sin^ with % G [0, vr/2), and 5, 021, and 031 represent the Dirac 
CP violation phase, and two Majorana CP violation phases, respectively. The CP phases 
vary from to 27r. 

Lastly let us derive a relation between Umns and U^. We define a unitary matrix Uh 
which diagonalizes DhU]^U'^D]^jUj^UjiDh as 



^ mi ^ 



ma 
mg 



(15) 



Then Umns is related to Ul and Uh as 

Umns = UlUh. (16) 

Note that, while Uh depends on the mass and mixing of the right-handed neutrinos, it is 
independent of the mixing of the lepton doublets, U^. This relation is important when 
we consider leptogenesis. 

B. Random matrix and measure 

The neutrino mass anarchy assumes the basis-independent random distribution of the 



matrices h and M [8|. Here we quote some of the results in Ref. jsl without derivation. 

Let us start with how to obtain a basis- independent random matrix, (/i)jj, with each 
element of order unity. We may generate a random number z for each element, uniformly 
distributed in a region of —1 < Re [2;] < 1 and — 1 < \m[z\ < 1. However, thus generated 
matrix is not basis- independent, as it changes its form under the ?7(3) rotation of the 
generations. In order to obtain a basis-independent random matrix, we need to impose 
Tr[/i/i^] < 1, which makes the distribution invariant under ?7(3). Similarly we can generate 
a random symmetric matrix X. 

Next question is the distribution of the eigenvalues /ii,/i2; and /is, and the mixing 
matrices Ul and f/^j, which should be invariant under U(3). The invariant measure of h 
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is given by 

3 

dh = FD{h,h2,h) f\dhi^^^^ (17) 

with 

FD{h,MM) = {hl-hl)\hl-hl)\hl-hlfh,h2h,, (18) 

where dUi and dUR represent the Haar measure of Ul and Ur, respectively. The dip in 
the denominator mods out the three redundant phases; we can see this by noting that 
the decomposition ( ITU]) is not unique, and it is invariant under multiphcation oi Ul and 
Ur by a diagonal unitary matrix. Similarly, the measure of M is given by 

3 

= Fm(Mi,M2,M3) JJdM,df/jv (19) 
1=1 

with 

Fm(Mi,M2,M3) = {Ml - Ml){Ml - Ml){Ml - MDM^M^M^, (20) 

where dU]<q represents the Haar measure of U^. 

Note that, while the distributions of the mixing angles are determined uniquely by the 
U(3) invariance, the measure of the eigenvalues can, in general, depend on an additional 
factor that is invariant under ?7(3), such as Tr[/i/i^] or det[/i]. Throughout this paper 
we assume that there is no such additional factor in the measure; this is called "linear 
measure" in the literature. Later we will briefly comment on how our results may change 
if other measure is adopted. 

Lastly let us give the measure of the neutrino mass matrix in the see-saw mechanism. 
The linear measure is given by 

3 3 ^JJ rITT 

dhdM cx FD{h,h2,h)FM{Mi,M2,M3) l[dh,l[dMj ^ (21) 

1=1 j=i 

where dU^R is the Haar measure of Unr = U^jUr. The neutrino mass eigenvalues can be 
obtained by diagonalizing rriy = UlDhUj^j^DjjUNRDhU'lv'^. We may replace dllL with 
dUuNS by using f lT6|) . Thus, we can see that the mixing angles obey the invariant Haar 
measure of Umns, since there is no way to distinguish three generations. It is given by 

dUuNS oc dsl2dc'l^dsl^dS da2ida3i. (22) 
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From f|T6l) we can see that the Haar distribution of Umns arises from the U(3)-invariance 
of Ul- Thus, even if the distribution of Uh is significantly distorted by e.g. leptogenesis, 
the Haar distribution of the Umns matrix remains intact. 

For practical purposes, it is useful to parametrize the light neutrino mass matrix as 

= C^D^^^Gv^ (23) 

where Mi obeys the linear measure distribution Fm{Mi, M2, M^) and G is a complex 
valued 3x3 random matrix generated as explained above. In this method, the neutrino 
Yukawa couplings of the right-handed neutrino mass eigenstates are directly obtained as 
the matrix G, and so, it is convenient when we discuss leptogenesis. 

Since the light neutrino mass matrix is obtained by the product of several random 
matrices, its mass eigenvalues exhibit a mild hierarchy. In particular, the ratio of mass 
squared difference, R = {m\ — m\) / {m\ — ml), can be naturally as small as the observed 
value ~ 1/30 for the normal hierarchy. See Fig. [1] for the distributions of the mass 
squared differences (Am|2,Am2i), -R, (mi,m2,m3), mi (to be defined in Eq. fl28|) ). the 
baryon-to-photon ratio r^^, and mgg (to be defined in Eq. (132|) ). Here we have assumed 
that the typical values of hi^ and Xij are of order unity, and imposed the constraints 
Tr[/i/it] < 1 and Tr[XXt] < 1, and fixed Mq to be 10^^ GeV. One can see from the 
figure that the distributions are consistent with the observation. In particular, the mild 
hierarchy R ~ 1/30 is nicely explained. Note that the distribution of R remains intact 
even if we change the typical values of the neutrino Yukawa matrix and Mq, while the 
neutrino mass distribution is affected. This is no longer the case when leptogenesis is 
taken into account, and both R and mj sensitively depend on the reheating temperature. 

We can see from the figure that the normal hierarchy is preferred since Amgg tends to 
be larger than Am^i. In our notation, the inverted hierarchy is realized if Am32 ^ Am^i, 
and in order to be consistent with the observations, the distribution of Am32 and Am^i 
should overlap with the left (blue) and right (red) dashed vertical lines, respectively. We 
can see that, although rare, the inverted hierarchy is indeed possible in the neutrino mass 
anarchy hypothesis. 
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FIG. 1: The distributions of the mass squared differences (Am32 (red), Am^]^ (blue)), their 
ratio R, the neutrino masses (mi (blue),m2 (green), ma (red)), mi (to be defined in Eq. (f28l) ). 
the baryon-to-photon ratio r]B, and niee (to be defined in Eq. (I32p ). The vertical lines represent 
the observed mass squared differences R ^ 0.03 (normal hierarchy) and R ^ 30 (inverted 
hierarchy) . 



We will study how the distributions are affected if we impose the leptogenesis in the 
next section. 
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III. LEPTOGENESIS AND ANARCHY 



We are interested in the conditional distribution of the neutrino parameters where 
leptogenesis works successfully. In the leptogenesis scenario, the lepton asymmetry is 
generated by the out-of-equilibrium decay of right-handed neutrinos. In the thermal 
leptogenesis scenario with zero initial abundance, the right-handed neutrino is generated 
thermally by inverse decay and scattering processes in thermal plasma. On the other 
hand, the right-handed neutrino is generated non-thermally by the inflaton decay in the 



non-thermal leptogenesis scenario 20|, |21| . We will see that this distinction is crucial in 



the neutrino mass anarchy. To simplify our analysis, we focus on a case in which the final 



lepton asymmetry is predominantly generated by the decay of the lightest ri 



it-handed 



neutrino, Ni, and we do not consider effects of the flavored leptogenesis 



22 



23|.^ The 



resonant leptogenesis [24[ is disfavored in the neutrino mass anarchy, because the measure 
( fT9l) forces the right-handed neutrino masses to be apart from each other. 

As is well known, the reheating temperature after inflation is bounded below, Tr > 
10^(10^) GeV, in order for (non-)thermal leptogenesis to account for the observed baron 



asymmetry Ref. 25|. Thus the reheating temperature Tr is an important input parameter 
for leptogenesis, but its precise value is poorly known, and so, we treat Tr as a free 
parameter and see how the distribution of the neutrino parameters changes as we vary 
Tr. We do not take account of the prior distributions of Tr and the resultant baryon 
asymmetry, because they are likely distorted by the anthropic conditions if leptogenesis is 
responsible for the origin of matter. On the other hand, if both the neutrino mass anarchy 
and leptogenesis are realized in nature, the observed neutrino mass squared differences 
and the mixing angles should be typical in the conditional distribution (as long as the 
light neutrino masses are irrelevant to the existence of life). 



^ We adopt a basis in which the right-handed neutrino mass matrix is diagonahzed in this section. See 
discussion below for how we generate random matrices. 
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A. Preliminaries 



The decay rate of A'^i at tree level is given by 



ri(iVi ^H + Q = f i(iVi ^H^ + il) 



(24) 



ignoring the masses of the final states. The CP asymmetry ei of the decay of A^i reads 



26 



with 



/(x) = x(l + x^) log (1 + X ^)-x + 



X 



x^ — 1 



Then, the final baryon-to-photon ratio ?7_b is given by 



r]B ~ 3 eiK, 



(25) 



(26) 



(27) 



where agph represents the sphaleron conversion factor, and it is equal to 29/78 in the SM, 
~ 100 counts the relativistic degrees of freedom at the A^i decay, and k denotes the 
efficiency factor 25|]. We take account of both AL = 1 and AL = 2 wash-out processes 
in our analysis. For later use, let us also define the following parameter: 



mi = 



(/i/it) 



11 



Ml 



(28) 



which is proportional to the ratio of the decay rate of A'^i to the Hubble parameter at the 
temperature T = Mi. 



B. Invariant Haar measure of Umns 

Here we show that the distribution of Uuns is orthogonal to the parameters relevant 
for leptogenesis. Since the lepton asymmetry is generated by the decay of A'^i, it is Mi and 
hia that are especially relevant for leptogenesis. As we have seen, the neutrino Yukawa 
matrix {h)ia always appears in a form of {hh^) in the decay rate, the CP asymmetry, and 
the efficiency parameter. Since (hh^) is invariant under the U(3) rotation of ia (see ([8])), 
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FIG. 2: The solid lines represent the Haar measure distribution of the mixing angles and CP 
phases given by (j22p . The red dots represent the distributions when we impose a successful 
leptog enesis constraint: Mi max 

= 10^^ GeV and 10 ^'^ < tjb < W ^ . 

Ul is not constrained by the leptogenesis. Whatever constraints on Mi and/or [hh"^) are 
imposed, it does not affect the distribution of Ul. (The distribution of Uh is affected.) 
Considering the relation between U^f^s and Ul given by f|T6|) and the translational invari- 
ance of the Haar distribution, we conclude that the distributions of the neutrino mixing 
angles as well as Dirac and Majorana CP violation phases are independent of leptogenesis. 
Thus, the measure of Umns is still given by (!22l) even if we require successful leptogene- 
sis. We have confirmed this property numerically as well. In Fig. [21 one can see that the 
mixing angles and the CP violation phases are determined by the Haar measure of U(3). 
This is a good news: the success of the neutrino mass anarchy about the mixing angles is 
maintained, even if we require the successful leptogenesis. 

C. Thermal leptogenesis 

In thermal leptogenesis with zero initial abundance, A^i is thermally produced by the 
inverse decay and scattering processes, and its out-of-equilibrium decay generates lepton 
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asymmetry, which is converted to the baryon asymmetry. Requiring successful thermal 
leptogenesis affects the original distribution of the neutrino parameters in two ways. First 
we consider the effect of the reheating temperature Tr. To simplify our argument, we 
assume that the radiation dominated Universe started with temperature Tr, and there 
was no thermal plasma before the reheating. In the case of the usual exponential decay, 
there is thermal plasma even before the reheating. However, even if the right-handed 
neutrinos are produced before reheating, they will be diluted by the subsequent entropy 
production. Such a crude approximation is sufficient for our purpose. 

If Tr is close to Mo(= 10^^ GeV), all the three right-handed neutrinos are thermalized, 
and the leptogenesis is possible. The distribution of the baryon asymmetry is shown 
in Fig. [U The conditional distribution is obtained if we limit ourselves to the neutrino 
parameters leading to the correct amount of the baryon asymmetry. 

If Tr <^ Mq, on the other hand, the three right-handed neutrinos tend to be too heavy 
to be produced. So, for most of the neutrino parameters, thermal leptogenesis does not 
work. However, although rare, the lightest right-handed neutrino can be light enough, 
by chance, to be thermally produced. So, successful leptogenesis is possible only in such 
subset Si satisfying Mi < -zTr: 

5i: Ml <2;Tfi<M2<M3, (29) 

where M2 and M3 are comparable to Mq. In the weak washout regime, z ^ 1, while z is 



about 4 — 6 in the strong washout regime 



25|. In Si, the distribution of Mi is peaked at 



Ml ~ zTji. It is clear that the neutrino mass distribution in Si is far from the observed 
one. That is to say, one of the light neutrino will be much heavier than the sub-eV scale. 
This can be understood by noting that the right-handed neutrino mass appears in the 
denominator in the see-saw formula for the light neutrino mass ( IT2|) . Therefore, as long 
as the neutrino Yukawa couplings are of order unity, one of the light neutrino masses will 
become much heavier than the other two. 

In order to avoid this problem, the neutrino Yukawa couplings, hia, must be suppressed. 
Such suppression of the Mi and hia can be easily realized by a simple U(l) flavor symmetry 
under which Ni is charged. However, as we are considering the neutrino mass anarchy 
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without any flavor symmetries in the neutrino sector, we need some other explanation. 
Intriguingly, such suppression is actually required in thermal leptogenesis. This is because 
of the wash-out effect: the neutrino Yukawa coupling hia must be suppressed because 
otherwise the resultant lepton asymmetry would be erased efficiently. Thus, the successful 
leptogenesis selects the following subset: 

^2 : < 1. (30) 

a 

This is the second effect of thermal leptogenesis. The analytic and numerical estimate 
of the upper bound can be found in e.g. Ref. {25]. We will see that the upper bound 
can be expressed in terms of fhi as rhi < 0(0. 1) eV. (Note that we have not imposed the 
observed mass squared difference.) 

To summarize, successful thermal leptogenesis selects the subset Sif]S2, where the 
hierarchical mass spectrum for the right-handed neutrinos. Mi ^ M2, M3, as well as the 
suppressed neutrino Yukawa couplings, \hia\ <^ 1, are realized. Such feature is genuinely 
emergent from the anarchy and thermal leptogenesis: no flavor symmetry is required. 

Before proceeding, let us describe our strategy. If the neutrino mass anarchy and the 
leptogenesis are indeed realized in nature, the observed neutrino parameters ([1]) and ([2]) 
should be typical ones in the conditional distribution. Since the mixing angles obey the 
Haar distribution, we will focus on the distribution of the neutrino masses. To this end, 
we will study the conditional distribution of the neutrino masses for which thermal lep- 
togenesis works successfully with given Tr. Note that we do not impose the observed 
values ([2]). If the observed mass squared differences are not typical in the obtained dis- 
tribution, we conclude that such a framework is disfavored; either the true reheating 
temperature should be different, or the assumptions of the neutrino mass anarchy and/or 
thermal leptogenesis are wrong. In this paper, we do not estimate the goodness of flt, 
since the qualitative understanding is fully adequate for our purposes. We leave a detailed 
statistical test for future work. 

Now let us go into details and see how the distributions change as we vary Tr. In 
Fig.[T]we have shown the original distribution. Now we impose the successful leptogenesis; 
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namely, we require the baryon asymmetry to be in the following range: 



.4 



5 X 10"'" < r]B <7 X 10 



-10 



(31) 



In Fig. |3] we show the distribution of the mass squared differences (Am32,Am2i), R, 
(mi,m2,m3), and rhi. We have generated one million sets of random matrices satisfying 
(l3Tll . Note that here we do not impose any constraint on Mi. This corresponds to the 
case of a high reheating temperature, Tr ~ Mq. We can see that the typical value of R 
decreases by one order of magnitude, compared to the original distribution in Fig. [1] This 
can be understood as follows. Since the original distribution of the baryon asymmetry 
is peaked around rjs ~ 10~^ (see Fig. [1]), we need to suppress the baryon asymmetry to 
satisfy Eq. (l3Tll . This can be achieved by decreasing the value of Mi, which leads to the 
increase of rhi. (Note that the value of rhi is bounded above, rhi < (9(0.1) eV, for the 
successful leptogenesis.) As a result, the mass of the heaviest neutrino, ms, tends to be 
heavier, suppressing R. We can see that the distribution of Am^g is in tension with the 
observations. 

Now let us consider a case of Tr ^ Mq. In this case, there is an upper bound on Mi 
for A'^i to be thermally produced. See Eq. ( I29l) . To simplify our analysis, we introduce a 
cut-off on Ml, Ml < Mi^rnax, whcrc the maximum value Mi^^ax is comparable to Tr.^ 

We show in Fig. H] the distributions of the neutrino parameters for Mi^^ax = 10^^ GeV. 
We can see that the upper bound on rhi is saturated at about O(0.1)eV, which results 
in the relatively heavy > 0.1 eV. As a result, the distribution of is peaked at 

~ 0.1 eV^ in strong tension with the observations. Also, the distribution of R is peaked 
below ~ 10~^ in slight tension with the observations. We emphasize here that this tension 
cannot be removed by simply changing the typical scale of Mq. If we increase Mq, the 
distribution of Am2i decreases almost in proportion to I/Mq, because the effect of m2 



^ The reason why we do not impose the observed value, rjB — (6.19 ± 0.15) x 10^^" [3], is to increase 
the number of random matrices satisfying the above criterion. Our main purpose here is to obtain the 
qualitative understanding of how the distributions of the neutrino parameters are modified by imposing 
successful leptogenesis. 

^ This is an approximation because the coefficient z in (|29p depends on the neutrino Yukawa couplings. 
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Logio(m/eV) Logio(rni/eV) 



FIG. 3: Same as Fig.[Tl except that we have required successful leptogenesis (f3T]l . No constraint 
on Ml is imposed. 



and nil on the leptogenesis is mild. On the other hand, the distribution of ms and 
therefore of Am^j does not change significantly, because it is determined by the balance 
between leptogenesis and random matrices of order unity. Thus, the distribution of R goes 
toward even smaller values, and the tension actually gets severer. For instance, if we take 
Mo = 3 X 10^^ GeV, the distribution of Am^^ can be in agreement with the observation, 
while that of R is peaked at about 10~^. If we decrease Mq, on the other hand, Am^i 
increases while Amgg remains almost the same. As a result, the mass squared difference 
Am2i will be much larger than the observed values. For instance, both Am^i and Am|2 
are about two orders of magnitude larger than the observed ones for Mq = 3 x 10^'^ GeV. 

The situation changes when we consider Mi^^ax ^ 0(10^^) GeV, for which it becomes 
difficult to generate the right amount of the baryon asymmetry and the upper bound 
on fhi from the wash-out effect becomes tighter, fhi < 0.1 eV. The distributions of the 
neutrino parameters are shown in Fig. |5]for Mi^max = 5 x 10^° GeV. One can see that the 
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Logio(m/eV) Logio(mi/eV) 

FIG. 4: Same as Fig. [3l except that we have imposed a constraint Mi < 10^'^ GeV. 

distribution of R is peaked at about 1/30, and also the mass squared differences agree 
very well with the observations. We also show the distribution of rriee- We have confirmed 
that the situation is similar for Mi max = 10^^ GeV. 

The results for Mi ^ax = 10^*^ GeV are shown in Fig. [6l One can see that the distribu- 
tions are consistent with the observations, although the constraint on rhi becomes even 
tighter and the distribution of R becomes broader. 

To summarize, the neutrino mass distribution nicely explains the observed neu- 
trino mass squared differences for Mi^max ~ (9(10^°^^^) GeV, or equivalently, Tr = 
0(1O^~^°) GeV in terms of the reheating temperature. On the other hand, the neutrino 
mass distribution is in tension with the observations for Tr higher than (9(10^^) GeV, 
and in particular the tension is significant for Tr around (9(10^^"^'^) GeV. The solution 
of Tr = (9(10^) GeV is particularly interesting from the point of view of inflation model 
building, dark matter, and the recent indications of a 125 GeV SM-like Higgs boson, which 
will be discussed in SeclIVi 
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We close this subsection by briefly discussing the distribution of ruee defined by 

3 



rriee = ^{UMNsYeim, 



(32) 



■t=l 



to which the amplitude for neutrinoless double beta decay (Oz//3/3) is proportional. The 



best limit for Ge is lmee| < 0.35 eV |29| , and the recent constraints for Xe are given 
by KamLAND-Zen [Soj and EXO-200 [aij experiments as, \m^e\ < (0.3 - 0.6) eV and 



< 



(0.14 - 0.38) eV at 90%C.L., respectively. 



20 




FIG. 6: Same as Fig. [3l except that we have imposed a constraint Mi < 10^° GeV. 



From the observed values of the mixing angles ([T]), one finds 

rriee ^ 0.67mi + CSOmse^^^^ + ^mm^e''^^^ , (33) 

implying that mee is very small for the neutrino masses with normal hierarchy. Indeed, 
iriee typically lies in the range of a few 10"'^ eV in the neutrino mass anarchy, and thus is 
below the reach of current experiments. Here we have used the fact that the distribution 
of Umns is subject to the invariant Haar measure and is independent of leptogenesis, 
and that the neutrino mass eigenvalues obey the conditional distribution where thermal 
leptogenesis works successfully and the observed mass squared differences are realized. 



D. Non-thermal leptogenesis 

It is not known how the reheating proceeds, because the coupling of the infiaton with 



the SM particles are poorly constrained. In a class of infiation models, the infiaton 



3^ or 
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waterfall field [33| is identified with the U(1)b-l Higgs boson, which is naturally coupled 
to the right-handed neutrinos to generate a large Majorana mass. Then the right-handed 
neutrinos are produced by the inflaton decay, and non-thermal leptogenesis takes place 



if the reheating temperature is lower than Mi 20|, l2l(| . The right amount of the baryon 
asymmetry can be created at a low reheating temperature, Tr > 10^ GeV, since the 
wash-out effect is suppressed. 

Let us consider non-thermal leptogenesis with the neutrino mass anarchy hypothesis. 
Suppose that the inflaton mass is much smaller than Mq. Then, the typical mass 
spectrum will be 

TK<Mi<m^< M2<Ms. (34) 

However, since the wash-out effect is weak at such low Tr, the constraint on the neutrino 
Yukawa coupling is much weaker than thermal leptogenesis. Therefore, the contribution 
of A^^i to the light neutrino mass will be signiflcantly larger than those of N2 and N3, and 
the resultant mass spectrum is 

mi <m2 <^ rris, (35) 

leading to an unacceptably small value of R, in contradiction with the observations. 
This problem can be avoided if the inflaton mass is heavier than or comparable to Mq: 

>Mo. (36) 

For the reference value of Mq = 10^^ GeV, this inequality is met only for a limited class 
of inflation models such as a smooth hybrid inflation js^. Note that this problem can be 
avoided in a flavor model in which the lightest right-handed neutrino is charged under a 
U(l) flavor symmetry. 



IV. DISCUSSION AND CONCLUSIONS 

We have so far flxed Mq = 10^^ GeV. Let us briefly discuss what happens if other 
values of Mq are chosen. First note that normal mass hierarchy is preferred compared to 
inverted hierarchy or degenerate spectrum in the neutrino mass anarchy. Thus, the mass 
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squared differences (Am32,Am2i) are approximately given by {ml^m^ for most of the 
cases, and tliese two mass scales must agree well with the observations. For Mq many 
orders of magnitude smaller than 10^^ GeV, it is therefore necessary to change the typical 
value of the neutrino Yukawa couplings in order to explain the observed neutrino mass 
squared differences. In particular, this is necessary because otherwise m2 (and therefore 
Am2i) would be in tension with the observations. For example, if Mq is suppressed by 
a factor of 10^, we should impose Tr[/i/it] < 10-2 instead of Tr[/i/it] < 1. On the other 
hand, Mq cannot be much larger than 10^^ GeV because the neutrino Yukawa couplings 
are required to be too large for perturbative computations. In these cases it is natural 
to introduce a common flavor charge of the right-handed neutrinos to explain the typical 
values of Mq and hi^. This is an interesting possibility, but it requires computationally 
expensive parameter scan. So, here let us focus on the neutrino Yukawa couplings of order 
unity and briefly mention how the distribution changes for a slightly different value of 
Mq. We have studied different values of Mq around 10^^ GeV, and found that the allowed 
region of Mi^^nax decreases as Mq deviates from 10^^ GeV. For instance, the distributions 
are in agreement with the observations for Mi^max ~ 3 x 10^^ GeV — 10^^ GeV, if we 
take Mq = 3 x lO""^^ GeV. As we increase Mq further, Am2i tends to be too small to 
account for the observed value. On the other hand, if we take Mq = 3 x 10^^ GeV, 
Am2i tends to be large, and Mi ~ 10^° GeV gives distributions barely consistent 
with the observations. Thus, for Mq = 3 x 10^^ GeV ~ 3 x 10^^ GeV, the neutrino mass 
distributions agree well with the observations for Mi = 10^° GeV ~ 10^^ GeV, or 
equivalently, Tr = 0(10^-") GeV. 

The SM with three right-handed neutrinos has been considered in our analysis, but 
it is straightforward to extend it to the supersymmetric (SUSY) framework. Our main 
conclusion still holds in this case. Interestingly, the ATLAS and CMS collaborations have 
recently pro vided hints for the existence of a SM-like Higgs particle with mass about 



125 GeV 



35| . The relatively light Higgs boson mass suggests the presence of new physics 



at scales below the Planck scale 



. In SUSY extensions of the SM, a 125 GeV Higgs 



mass can be explained without invoking large stop mixing if the typical sparticle mass is at 
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0(10) TeV or heavier. Among various possibilities of the SUSY breaking mediation mech 
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anisms, the simplest one is the anomaly mediation with a generic Kahler potential 
Then the Wino is likely the lightest SUSY particle (LSP), and therefore a candidate for 
dark matter if the R-parity is conserved. For the gravitino mass of O(100)TeV, the 
thermal relic abundance of the Wino is too small to account for the observed dark mat- 
ter density. The correct abundance can be naturally realized by the gravitino decay, if 
Tr = (9(10^) GeV. Interestingly enough, with this reheating temperature, the thermal 
leptogenesis is possible. Furthermore, we have seen in the previous section that the neu- 
trino mass spectrum is typical in the conditional distribution at Tr = O{10^^^^) GeV for 
Mq ^ 10^^ GeV when the neutrino mass anarchy is assumed. Thus, the 125 GeV SM-like 
Higgs boson, the LSP dark matter produced by the gravitino decay, thermal leptogene- 
sis, and the neutrino mass anarchy point to Tr = 0(10^"^'') GeV. Such a coincidence is 
interesting and even suggestive. 

We have assumed that there is no flavor symmetry which distinguishes three genera- 
tions of the neutrinos. The observed hierarchical spectrum of quarks and charge leptons, 
on the other hand, can be nicely explained by the flavor symmetry under which only 
10-plets are charged while 5*-plets are neutral in the language of the SU(5) Irre- 
spective of whether the flavor symmetry is a true symmetry or an emergent one, this 
is consistent with the SU(5) GUT. We also note that some of the problems outlined in 
the previous section (e.g. the difficulty in the non-thermal leptogenesis) can be easily 
solved if there is a flavor symmetry under which the right-handed neutrinos are charged. 
It should be emphasized that such flavor symmetry does not affect the see-saw formula 
for the light neutrino mass. In fact, it was shown that the see-saw mechanism is robust 
against splitting the right-handed neutrino masses in this way 



One of the important assumptions in our analysis is that both hia and Xij obey random 
distribution of order unity. The typical value can be different from order unity by assigning 
a common flavor charge on three generations of A^, and/or i^, or by an extra dimensional 
set-up. In this sense the neutrino mass anarchy and the conventional flavor symmetry are 
compatible. See Ref. j40j] for the recent study of the neutrino mass anarchy with a certain 
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flavor symmetry [41]. 

In our analysis we have assumed the linear measure of hia and Xij. It is in principle 
possible to adopt another measure which depends on a U(3)-invariant factor, such as 
Tr[/;./i^] or det[/i]. In this case, while the mixing angles and CP violation phases are still 
determined by the invariant Haar measure, the distributions of the neutrino masses are 
generically modified. However, we believe that our results are robust against changing 
the measure to some extent. This is because the peak position of is determined by 
the balance between the randomness and the wash-out effect. To see this, let us consider 
thermal leptogenesis with a reheating temperature lower than the typical right-handed 
neutrino mass. Recall that the successful thermal leptogenesis selects the subset 5i f]52, 
where Mi <^ M2,M3 and \hia\ <^ 1 are realized. Therefore, for given Tr, the peak of 
the distribution of rris is not sensitive to the details of the adopted measure, as long as 
the conditional distribution of Mi is peaked near and the wash-out bound on the 
\hia\ is saturated. Suppose, for instance, that the measure is an increasing function of 
the eigenvalues of h and M; then the upper bounds on Mi and hia are considered to be 
saturated because it requires severe fine-tunings to realize smaller values. Therefore, as 
long as the measure of h^a and Xij satisfy such property, the distribution of ms is not 
modified significantly. The distribution of mi and m2 are not sensitive to the leptogenesis, 
and so, our results are valid if the measure is such that the predicted neutrino mass 
spectrum agrees with the observation when successful leptogenesis is not imposed. Of 
course, our results do not hold, for instance, if the measure favors a lighter right-handed 
neutrinos or smaller neutrino Yukawa couplings. 

The observed neutrino mixing angles and mass squared difference support both the 
neutrino mass anarchy and the see-saw mechanism. This suggests that the neutrino 
masses and mixings are irrelevant to the existence of life and therefore keeps its original 
distribution in the landscape. On the other hand, however, if the leptogenesis is respon- 
sible for the origin of matter, it selects a certain subset of the neutrino parameters, which 
may significantly distort the original distribution spoiling the success of neutrino mass 
anarchy. In this paper we have studied if the successful leptogenesis is possible together 
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with the neutrino mass anarchy hypothesis. We have found that the distributions of the 
neutrino mixing angles as well as Dirac and Majorana CP violation phases are determined 
by the invariant Haar measure of U(3), even if we impose successful Icptogcnesis. On the 
other hand, the neutrino mass spectrum is generically affected by leptogenesis. In the 
case of thermal leptogenesis, the mass spectrum for the right-handed neutrinos and the 
neutrino Yukawa matrix exhibit a certain pattern, as a result of the competition between 
random matrices with elements of order unity and the wash-out effect. The hierarchical 
mass spectrum for the right-handed neutrinos as well as suppressed neutrino Yukawa cou- 
plings hia are similar to that obtained by an approximate flavor symmetry. In a sense, 
the flavor symmetry is emergent. However, as we look into details, we have seen that 
the neutrino mass spectrum depends sensitively on the reheating temperature. The light 
neutrino mass distribution is consistent with observation if the reheating temperature is 
0(109-11) GeV for Mq = 3 x IQi^GeV ~ 3 x IQi^GeV. In particular, the solution of 
Tr — 0(10^) GeV is interesting in connection with the neutrahno dark matter produced 
by the gravitino decay and the 125 GeV SM-like Higgs boson. On the other hand, non- 
thermal leptogenesis is consistent with the observations only if the inflaton mass is heavier 
than the typical right-handed neutrino mass scale. Such a heavy inflaton mass can be 
realized in e.g. the smooth hybrid inflation. 
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